Genome-wide association studies for heat stress response in Bos taurus × Bos indicus crossbred cattle.

Heat stress is an important issue in the global dairy industry. In tropical areas, an alternative to overcome heat stress is the use of crossbred animals or synthetic breeds, such as the Girolando. In this study, we performed a genome-wide association study (GWAS) and post-GWAS analyses for heat stress in an experimental Gir × Holstein F2 population. Rectal temperature (RT) was measured in heat-stressed F2 animals, and the variation between 2 consecutive RT measurements (ΔRT) was used as the dependent variable. Illumina BovineSNP50v1 BeadChip (Illumina Inc., San Diego, CA) and single-SNP approach were used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene-transcription factor (TF) networks, generated from enriched TF. The breed origin of marker alleles in the F2 population was assigned using the breed of origin of alleles (BOA) approach. Heritability and repeatability estimates (± standard error) for ΔRT were 0.13 ± 0.08 and 0.29 ± 0.06, respectively. Association analysis revealed 6 SNP significantly associated with ΔRT. Genes involved with biological processes in response to heat stress effects (LIF, OSM, TXNRD2, and DGCR8) were identified as putative candidate genes. After performing the BOA approach, the 10% of F2 animals with the lowest breeding values for ΔRT were classified as low-ΔRT, and the 10% with the highest breeding values for ΔRT were classified as high-ΔRT. On average, 49.4% of low-ΔRT animals had 2 alleles from the Holstein breed (HH), and 39% had both alleles from the Gir breed (GG). In high-ΔRT animals, the average proportion of animals for HH and GG were 1.4 and 50.2%, respectively. This study allowed the identification of candidate genes for ΔRT in Gir × Holstein crossbred animals. According to the BOA approach, Holstein breed alleles could be associated with better response to heat stress effects, which could be explained by the fact that Holstein animals are more affected by heat stress than Gir animals and thus require a genetic architecture to defend the body from the deleterious effects of heat stress. Future studies can provide further knowledge to uncover the genetic architecture underlying heat stress in crossbred cattle.

Genome-wide association studies for tick resistance in Bos taurus × Bos indicus crossbred cattle: A deeper look into this intricate mechanism.

Rhipicephalus (Boophilus) microplus is the main cattle ectoparasite in tropical areas. Gir × Holstein crossbred cows are well adapted to different production systems in Brazil. In this context, we performed genome-wide association study (GWAS) and post-GWAS analyses for R. microplus resistance in an experimental Gir × Holstein F2 population. Single nucleotide polymorphisms (SNP) identified in GWAS were used to build gene networks and to investigate the breed of origin for its alleles. Tick artificial infestations were performed during the dry and rainy seasons. Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA) and single-step BLUP procedure was used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene transcription factors networks, generated from enriched transcription factors, identified from the promoter sequences of selected gene sets. The genetic origin of marker alleles in the F2 population was assigned using the breed of origin of alleles approach. Heritability estimates for tick counts were 0.40 ± 0.11 in the rainy season and 0.54 ± 0.11 in the dry season. The top ten 0.5-Mbp windows with the highest percentage of genetic variance explained by SNP markers were found in chromosomes 10 and 23 for both the dry and rainy seasons. Gene network analyses allowed the identification of genes involved with biological processes relevant to immune system functions (TREM1, TREM2, and CD83). Gene-transcription factors network allowed the identification of genes involved with immune functions (MYO5A, TREML1, and PRSS16). In resistant animals, the average proportion of animals showing significant SNPs with paternal and maternal alleles originated from Gir breed was 44.8% whereas the proportion of animals with both paternal and maternal alleles originated from Holstein breed was 11.3%. Susceptible animals showing both paternal and maternal alleles originated from Holstein breed represented 44.6% on average, whereas both paternal and maternal alleles originated from Gir breed animals represented 9.3%. This study allowed us to identify candidate genes for tick resistance in Gir × Holstein crossbreds in both rainy and dry seasons. According to the origin of alleles analysis, we found that most animals classified as resistant showed 2 alleles from Gir breed, while the susceptible ones showed alleles from Holstein. Based on these results, the identified genes may be thoroughly investigated in additional experiments aiming to validate their effects on tick resistance phenotype in cattle.